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Dear  Dr.  Nachman, 

Below  is  a  summary  of  our  research  effort  supported  by  the  aforementioned  grant. 

(1).  Theoretical  model  and  computer  program  for  3-D  laser  transmission/backscattering  through 
cirrus  clouds 


We  have  developed  a  3-D  radiative  transfer  model  based  on  the  successive-order-of-scattering 
approach  for  the  computation  of  airborne  continuous-wave  laser  transmission  and  backscattering 
through  3-D  inhomogeneous  high-level  clouds  in  both  plane-parallel  and  spherical  geometries. 
This  model  takes  into  account  the  relative  positions  of  aircraft,  target,  and  clouds  in  3-D 
atmospheres  and  is  formulated  on  the  basis  of  the  theory  of  radiative  transfer.  Subsequently,  3-D 
laser  transmission/backscattering  computer  programs  have  been  developed  for  potential 
incorporation  in  desk-top  decision-aid  models  to  determine  the  attenuation  of  high-energy  laser 
beam  through  3-D  cirrus  clouds  for  realistic  scenarios. 

In  conjunction  with  the  3-D  high  energy  laser  transmission  and  backscattering,  we  developed  an 
innovative  satellite  remote  sensing  approach  to  map  3-D  cirrus  ice  water  content  and  men 
effective  ice  crystal  size  fields  from  satellite  data  using  a  combination  of  the  conventional 
retrieval  of  cloud  optical  depth  and  particle  size  in  a  horizontal  plan  and  a  parameterization  of  the 
vertical  cloud  profile  involving  temperature  from  sounding  and/or  analysis.  Our  approach  has 
been  validated  in  part  by  two  impressive  cirrus  clouds  that  occurred  in  the  vicinity  of  northern 
Oklahoma  on  18  April  1997  and  9  March  2000  obtained  from  DOE’s  Atmospheric  Radiation 
Program  in  which  the  ice  crystal  size  distributions  were  collected  independently  from  collocated 
and  coincident  aircraft  optical  probe  measurements  under  the  field-of-view  of  AVHRR/NOAA. 

We  investigated  the  sensitivity  of  3-D  laser  transmission  and  backscattering  on  variation  in  the 
cirrus  cloud  optical  depth,  particle  size,  and  cloud-base  height,  as  well  as  the  positions  of  aircraft 
and  target.  We  show  the  significance  of  3-D  cloud  geometry  and  inhomogeneity  and  Earth’s 
spherical  geometry  on  the  transmitted  and  backscattered  laser  powers.  Moreover,  we  demonstrate 
that  the  3-D  cloud  fields  derived  from  satellite  remote  sensing  can  be  employed  for  the  3-D  laser 
transmission  and  backscattering  model  for  tactical  application.  Some  pertinent  results  have  been 
presented  in  the  attached  paper  to  be  published  in  Applied  Optics  in  September,  2006  (see  Figs.  3, 
4,  and  5). 
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Finally,  substantial  differences  are  shown  between  laser  transmissions  through  the  present  3-D 
model  and  previously  developed  2-D  model  (on  a  plane)  on  the  basis  of  numerical  simulations  in 
terms  of  direct  transmission,  forward-scattered  transmission,  and  backscattering.  Thus,  the  high 
energy  laser  transmission  program  involving  3-D  aircraft-target  geometry  must  account  for  the 
3D  inhomogeneity  of  cloud  fields  in  order  to  realistically  evaluate  the  transmitted  and 
backscattered  powers  for  tactical  application.  We  also  find  that  the  effect  of  Earth’s  curvature 
becomes  particularly  important  for  long  path  transmission  involving  the  airborne  laser  and 
moving  target  engagement. 

(2).  Remote  Sensing  of  Thin  Cirrus/ Aerosols  to  Support  Airborne  Laser  Transmission  Program 

We  developed  a  new  methodology  to  simultaneously  retrieve  the  optical  depths  of  aerosols  and 
thin  cirrus  clouds  over  the  oceans  by  using  the  MODIS/NASA  data.  This  approach  takes 
advantage  of  the  fact  that  the  visible  or  a  near-IR  window  reflectance  from  cirrus  can  be 
characterized  by  its  correlation  to  the  reflectance  from  a  near  IR  band  located  in  a  strong  water 
vapor  absorption  region.  The  cirrus  contribution  to  the  total  reflectance  can  be  first  removed, 
followed  by  the  retrieval  of  the  aerosol  optical  depth  from  the  residual  reflectance  in  these  bands. 
Iterations  have  been  developed  to  obtain  the  cirrus  optical  depth  by  the  principle  of  minimization 
between  observed  and  calculated  reflectance  values.  Examination  of  the  data  gathered  from  the 
recent  CRYSTAL-FACE  and  CLAMS  field  campaigns  demonstrates  that  our  retrievals  are 
consistent  with  both  in-situ  and  ground  based  measurements.  The  preceding  results  have  been 
published  in  Roskovensky  et  al.  (2004). 
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Remote  sensing  of  three-dimensional  cirrus  clouds  from 
satellites:  application  to  continuous-wave  laser 
aQ  a  atmospheric  transmission  and  backscattering 


K.  N.  Liou,  Szu-Cheng  Ou,  Yoshihide  Takano,  and  Jeffrey  Cetola 


A  satellite  remote  sensing  methodology  has  been  developed  to  retrieve  3D  ice  water  content  (IWC)  and 
mean  effective  ice  crystal  size  of  cirrus  clouds  from  satellite  data  on  the  basis  of  a  combination  of  the 
conventional  retrieval  of  cloud  optical  depth  and  particle  size  in  a  horizontal  plane  and  a  parameteriza¬ 
tion  of  the  vertical  cloud  profile  involving  temperature  from  sounding  and/or  analysis.  The  inferred  3D 
cloud  fields  of  IWC  and  mean  effective  ice  crystal  size  associated  with  two  impressive  cirrus  clouds  that 
occurred  in  the  vicinity  of  northern  Oklahoma  on  18  April  1997  and  9  March  2000,  obtained  from  the 
Department  of  Energy’s  Atmospheric  Radiation  Measurement  Program,  have  been  validated  against  the 
ice  crystal  size  distributions  that  were  collected  independently  from  collocated  and  coincident  aircraft 
optical  probe  measurements.  The  3D  cloud  results  determined  from  satellite  data  have  been  applied  to 
the  simulation  of  cw  laser  energy  propagation,  and  we  show  the  significance  of  3D  cloud  geometry  and 
inhomogeneity  and  spherical  atmosphere  on  the  transmitted  and  backscattered  laser  powers.  Finally,  we 
demonstrate  that  the  3D  cloud  fields  derived  from  satellite  remote  sensing  can  be  used  for  the  3D  laser 
transmission  and  backscattering  model  for  tactical  application.  ©  2006  Optical  Society  of  America 
OCIS  codes:  280.0280,  120.0280,  140.0140,  010.3310. 


1 .  Introduction 

The  time  series  of  backscattering  coefficients  derived 
from  the  ground-based  lidar  and  cloud  radar  returns1-3 
have  shown  significant  vertical  variability  in  ice  water 
content  (IWC)  and  ice  crystal  size.  Vertical  inhomoge¬ 
neity  of  the  ice  crystal  size  and  IWC  has  also  been 
observed  in  the  balloonborne  replicator  sounding  ob¬ 
servations.4  The  vertical  inhomogeneity  of  ice  crystal 
size  can  alter  the  radiative  heating  profile  in  cirrus 
clouds  and  hence  is  an  important  component  in  the 
radiation- climate  feedback  study.5  It  can  also  affect 
the  propagation  of  active  radiation  sources  in  cloud 
layers.  However,  the  present  satellite  remote  sensing 
technique  can  only  retrieve  vertically  integrated  opti¬ 
cal  and  microphysical  properties  of  cirrus  clouds.6-12 
Cirrus  clouds  are  high  clouds  that  are  present  at  all 
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latitudes  and  in  all  seasons  with  a  global  average  of 
approximately  20%-30%  and  more  than  60%-70%  in 
the  tropics.13  Some  of  these  high  cirrus  clouds  are  op¬ 
tically  thin  with  optical  depths  less  than  0.5  (Refs.  14 
and  15)  and  often  subvisual.  In  this  paper,  we  have 
developed  a  parameterization  approach  to  derive  3D 
cloud  IWC  and  ice  crystal  size  distributions  based  on 
satellite  data  and  provided  a  number  of  independent 
validations  employing  concurrent  and  simultaneous 
in  situ  ice  microphysics  measurements  from  the  air¬ 
craft  platform  under  the  satellite  passage. 

In  conjunction  with  high-energy  laser  propagation 
in  the  atmosphere,  we  have  presented  ID  and  2D  laser 
transmission  models  in  plane-parallel  and  spherical 
geometries  associated  with  the  aircraft- cirrus-target 
system.16-18  Using  the  successive  order-of-scattering 
approach  for  radiative  transfer,  direct  and  diffuse 
transmissions  were  formulated  specifically  for  homo¬ 
geneous  high-level  clouds,  in  which  scattering  and  ab¬ 
sorption  associated  with  aerosols,  water  vapor,  and  air 
molecules  were  accounted  for.  In  the  2D  laser  trans¬ 
mission  model,18  it  was  assumed  that  cloud  properties 
only  vary  in  the  vertical  plane  ( x—z  plane),  thus 
neglecting  the  variation  of  the  cloud  properties  in  they 
direction.  Moreover,  the  3D  nature  of  the  laser  beam 
has  not  been  accounted  for  in  the  2D  formulation.  To 
accurately  simulate  the  3D  laser  transmission- 
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backscattering  in  high-level  cloud  fields,  a  3D  ap¬ 
proach  must  be  taken.  In  association  with  the 
mapping  of  3D  cloud  fields  from  satellite  remote 
sensing,  we  formulated  and  developed  a  3D  laser 
transmission- backscattering  model,  taking  into  ac¬ 
count  the  Earth’s  spherical  geometry  and  the  3D 
finite  and  inhomogeneous  structure  of  high-level 
cloud  fields  to  investigate  the  significance  of  the  3D 
nature  of  clouds  on  laser  energy  propagation  in  the 
atmosphere. 

This  paper  is  organized  as  follows:  Section  2  pre¬ 
sents  a  parameterization  approach  to  the  mapping  of 
3D  inhomogeneous  cirrus  cloud  fields  from  satellite 
data  and  application  of  the  3D  cirrus  cloud  remote 
sensing  method  to  Atmospheric  Radiation  Measure¬ 
ment  (ARM)  Cloud  Intensive  Operational  Period 
AQ:  b  (IOP)  data  with  an  independent  validation  of  its  ac¬ 
curacy  by  using  in  situ  aircraft  microphysical  mea¬ 
surements.  A  description  of  the  basic  formulation  and 
analysis  of  3D  laser  transmission  and  backscattering 
through  cirrus  clouds  follows  in  Section  3,  including 
an  assessment  of  the  importance  of  3D  cloud  effects 
on  high-energy  laser  transmission.  Finally,  conclud¬ 
ing  remarks  are  given  in  Section  4. 

2.  Satellite  Mapping  of  Three-Dimensional 
Inhomogeneous  Cirrus  Clouds 

A.  Parameterization  Approach 

A  number  of  approaches  have  been  developed  for 
the  inference  of  the  horizontal  distributions  of  cirrus 
cloud  optical  depth  and  ice  crystal  size  from  satellites. 
These  approaches  used  wavelengths  at  which  absorp¬ 
tion  by  water  vapor  and  other  gases  is  minimal  and  at 
which  scattering  and  absorption  by  cloud  particles 
exhibit  maximum  sensitivity.  Reflectance  correlation 
of  a  nonabsorbing  wavelength  (e.g.,  a  visible  wave¬ 
length)  and  an  absorbing  wavelength  (e.g.,  wave¬ 
lengths  at  1.6,  2.2,  or  3.7  pm)  has  been  established  by 
radiative  transfer  calculations  a  priori  in  terms  of 
lookup  tables  for  the  simultaneous  retrieval  of  optical 
depth  and  ice  crystal  size.7-10  The  cloud  optical  depth 
by  definition  refers  to  the  vertically  integrated  value, 
while  the  particle  size  retrieved  is  related  to  its  mean 
property  in  the  vertical.  Thus,  the  present  satellite 
retrieval  techniques  are  limited  to  the  horizontal 
mapping  of  the  cloud  optical  and  microphysical  prop¬ 
erties.  By  combining  the  horizontal  mapping  of  cloud 
parameters  from  satellites  with  the  retrieval  of  the 
vertical  distributions  of  IWC  and  ice  crystal  size  pa¬ 
rameters  from  the  reflectivity  and  Doppler  velocity  of 
a  millimeter-wave  cloud  radar  (mmCR),  Liou  et  al. 19 
showed  that  3D  distributions  of  IWC  and  ice  crystal 
mean  size  can  be  inferred.  This  novel  approach,  how¬ 
ever,  is  restricted  to  single  points  where  cloud  radar 
observations  are  available.  In  the  following,  we  de¬ 
velop  a  parameterization  approach  to  determine  the 
regional  3D  cloud  structure  entirely  from  satellite 
data  with  a  global  perspective.  We  will  also  present 
limited  validations  of  the  present  approach  based  on 
available  cloud  microphysics  measurements.  Our  ap¬ 
proach  focuses  on  global  application  and  the  use  of 


satellite  data  and  is  not  intended  to  account  for  local 
small-scale  cloud  variability  that  can  be  simulated 
from  stochastic  cloud  models.20’21 

We  first  let  the  horizontal  distribution  of  the  op¬ 
tical  depth  and  mean  effective  ice  crystal  size  re¬ 
trieved  from  satellite  data  be  denoted  by  t(x,  y)  and 
DE*fic,  y),  respectively.  We  define  DE  in  the  form 


DE  =  J  D2Ln(L)dL  j  j  DLn(L)AL,  (2.1) 

where  D  and  L  denote  the  width  and  the  maximum 
dimension  of  a  nonspherical  ice  crystal,  respectively, 
and  n(L)  is  the  size  distribution  in  terms  of  L.  The 
preceding  definition  of  the  mean  effective  size  is 
closely  related  to  the  ratio  of  the  total  ice  crystal 
volume  to  the  total  ice  crystal  cross-sectional  area 
and  is  applicable  to  irregular  and  complex  ice  crystal 
shapes.  Based  on  cloud  microphysics  observations 
and  radiative  transfer  parameterizations,  ice  crystal 
size  can  be  defined  in  terms  of  cloud  temperature, 
IWC,  and  optical  depth. 22-26  Moreover,  ice  crystal  size 
is  approximately  proportional  to  the  1/3  power  of 
IWC.22  Thus,  on  the  basis  of  dimensional  analysis,  we 
have  developed  the  following  parameterization  for 
the  3D  distribution  of  DE  in  the  form22 

DE(x,  y,  z )  =  c{j(x,  y)/[Az(ot  +  p/(DE(z)>) 

X  (IWC(z))]}1/3(DE(z)),  (2.2) 

where  t  is  the  visible  optical  depth,  c  is  a  proportional 
factor,  A z  is  the  cloud  thickness,  a  and  p  are  empirical 
constants,  and  the  terms 

(DE(z))  =  2  cn(Tc(z)  -  273)",  (2.3) 

n-0 

(IWC(z)>  =  exp(-7.6  +  4  exp{-0.2443  X  1CT3 
X  [253  -  Tc(z)]}2'445) 
for  Tc(z)  <  253  K.  (2.4) 

In  Eqs.  (2.3)  and  (2.4),  the  local  values  of  DE  and  IWC 
denoted  by  (•)  are  related  to  the  local  cloud  tempera¬ 
ture  Tc.  Information  about  Tc  (as  a  function  of  height) 
can  be  obtained  via  local  sounding.  Temperature  data 
may  also  be  obtained  from  the  conventional  retrieval 
using  vertical  sounding  channels27  or  from  the  Na¬ 
tional  Centers  for  Environmental  Prediction  (NCEP) 
global  analysis  product.28 

The  proportionality  factor  c  can  be  determined 
by  performing  a  vertical  averaging  of  Eq.  (2.2),  and 
equating  the  resulting  vertical  mean  DE(x,  y)  with 
the  satellite-retrieved  mean  effective  size,  DE:|:(x,  y), 
as  follows: 

c  =  DE*(x,  y)(T(z,  y)/{Az[a  +  0/DE*(:r,  y)] 
x  (IWC)})-1/3(DE>-1.  (2.5) 
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Finally,  the  3D  IWC  can  now  be  expressed  via  a 
radiative  transfer  parameterization  by29 

IWC(*,  y,  z)  =  t(x,  y)/{ Az[a  +  3/DE(x,  y,  z)]}. 

(2.6) 

The  methodology  for  the  remote  sensing  of  3D  cloud 
fields  outlined  above  is  based  on  a  number  of  physical 
approximations  and  local  microphysics  observations 
and  does  not  require  the  use  of  ground-based  lidar- 
radar  data.  Obviously,  its  success  requires  inde¬ 
pendent  validation.  We  present  two  case  studies  to 
illustrate  the  fundamental  principle  and  numerical 
procedure  of  the  present  approach.  For  these  two 
cases,  collocated  and  coincident  cirrus  cloud  micro¬ 
physical  data  available  from  in  situ  aircraft  measure¬ 
ments  during  the  same  time  period  of  the  satellite 
passage  were  used  for  an  independent  check  of  the 
retrieval  results.  It  should  be  noted  that  collocated  and 
coincident  satellite  and  in  situ  aircraft  data  requires  a 
carefully  planned  field  campaign  and  its  availability  is 
extremely  limited. 

B.  Validation  Using  Department  of  Energy-Atmospheric 
AQ:  e  Radiation  Measurement  Data 

To  cross-check  the  validity  of  the  parameterization 
approach  for  the  mapping  of  3D  cloud  fields  based 
on  satellite  data,  we  have  selected  Advanced  Very 
High-Resolution  Radiometer-National  Oceanic  and 
Atmospheric  Administration  (AVHRR-NOAA)  data 
gathered  at  the  ARM  Southern  Great  Plains  (SPG) 
Cloud  and  Radiation  Testbed  (CART)  site  in  the 
AQ:  f  northern  Oklahoma  area  for  single-layer  cirrus 
clouds  that  occurred  on  18  April  1997  and  9  March 
2000  during  two  Spring  Cloud  IOPs.  Retrieval  and 
analysis  were  carried  out  to  construct  3D  IWC  and 
DE  fields  from  the  AVHRR  data,  which  were  then 
compared  to  collocated  and  coincident  in  situ  mea¬ 
surements  of  the  ice  crystal  size  distribution  and  IWC 
independently  derived  from  the  optical  probes  on¬ 
board  the  University  of  North  Dakota  Citation. 

During  the  ARM  Spring  Cloud  IOP  at  the  SGP  on 
1-26  March  2000,  there  was  a  significant  stratiform 
cirrus  cloud  episode  occurring  on  9  March.  Well- 
defined  cirrus  clouds  started  to  appear  over  the  Cen¬ 
tral  Facility  site  after  1500  UTC,  reaching  peak 
activity  around  2100  UTC,  before  breaking  up.  Al¬ 
though  NASA/Terra  was  launched  three  months  be¬ 
fore  this  date,  Moderate  Resolution  Imaging  Spectro- 
aq-.g  radiometer  (MODIS)  data  for  this  date  was  not 
available  owing  to  insufficient  calibration  accuracy. 
However,  NOAA-14  passed  over  the  SGP  site  at  2206 
UTC  with  a  good  viewing  angle  and  AVHRR  data  was 
available  for  3D  cloud  studies.  The  horizontal  map¬ 
ping  of  cloud  parameters  is  based  on  the  correlation 
approach  for  the  0.63  |xm  reflectance  and  the  3.7  pm 
total  radiance  as  functions  of  reference  optical  depth 
and  DE  constructed  from  plane-parallel  radiative 
transfer  calculations  in  terms  of  lookup  tables  by  pre¬ 
scribing  the  single-scattering  properties  for  ice  par¬ 
ticles  a  priori ,9  It  should  be  noted  that  the  satellite 


retrieval  procedure  is  intricate  and  the  problem  of  3D 
inhomogeneity  in  radiative  transfer  is  one  of  many 
issues  in  developing  a  retrieval  algorithm.  The  dis¬ 
tribution  of  satellite  data  points  on  the  constructed 
charts  for  a  given  satellite-Sun  geometry  can  be  con¬ 
verted  to  optical  depth  and  DE  using  a  minimization 
method.  The  9  March  case  reveals  that  the  cirrus 
pixels  for  this  date  consisted  of  a  range  of  optical 
depths  from  0  to  6.  For  optical  depths  smaller  and 
larger  than  1,  the  DE  varies  between  30  and  130  pm 
and  between  70  and  130  pm,  respectively. 

A  3D  cloud  domain  containing  a  40  km  X  300  km 
rectangular  strip  was  chosen  at  which  the  center  co¬ 
incides  with  the  ARM  Central  Facility  site,  while  its 
orientation  is  along  the  wind  direction.  The  vertical 
profiles  of  temperature  and  wind  speed  and  direction 
from  soundings  launched  at  2030,  2206,  and  2329 
UTC  were  used  to  determine  the  vertical  IWC  and  ice 
crystal  size  using  the  cloud-level  temperature  sound¬ 
ing  and  following  the  cloud-top  height  determination 
developed  by  Hutchison  et  al.30  In  this  case,  the  peak 
wind  speed  at  the  cloud  level  is  approximately  40 
m/s,  and  the  wind  direction  is  relatively  constant  and 
maintained  at  approximately  250°,  approximately  in 
the  west-southwest  direction.  For  illustration  pur¬ 
poses,  the  3D  IWC  and  DE  results  are  presented  in 
xy,  yz,  and  xz  planes  over  a  120  km  X  40  km  X  9  km 
AVHRR  data  domain,  as  shown  in  Fig.  1  (left  panel). 
The  horizontal  and  vertical  resolutions  are  of  the 
order  of  1  and  0.1  km,  respectively.  Substantial  vari¬ 
ability  in  IWC  and  DE  is  shown  in  the  xz  plane,  as 
compared  to  those  in  the  yz  plane  associated  with  the 
prevailing  wind  direction.  Also  noted  is  the  increas¬ 
ing  DE  toward  the  cloud  base.  The  black  areas  indi¬ 
cate  that  AVHRR  data  was  not  available. 

The  satellite-inferred  3D  IWC  and  DE  fields  were 
compared  to  the  ice  crystal  size  distribution  and  IWC 
independently  derived  from  the  simultaneous  in  situ 
measurements  by  2D  optical  probes  on  board  the  Uni¬ 
versity  of  North  Dakota  Citation.  The  2D  probes  mea¬ 
sured  ice  crystal  sizes  in  25  p,m  increments  (type  C, 
from  25  to  100  p,m)  and  in  100  |xm  increments  (type 
P,  from  25  to  2000  (xm).31  In  this  case,  the  aircraft 
descends  from  9.5  to  6.5  km  and  the  flight  track 
between  2132  and  2223  UTC  consists  of  a  12-loop 
Lagrangian  spiral  descending  pattern  centered  at  the 
ARM  central  facility  site.  Each  loop  took  about  4  min 
to  complete,  and  the  aircraft  descended  0.25  km  per 
loop.  Sampling  of  the  ice  crystal  size  was  taken  every 
5  s.  Following  the  method  developed  in  Ou  et  al.,22  we 
computed  the  in  situ  IWC  and  DE  for  each  loop. 
Based  on  the  soundings  of  local  wind  speed  and  di¬ 
rection,  we  also  determined  the  latitude  and  longi¬ 
tude  of  the  satellite’s  pixels  that  correspond  to  the 
location  of  each  in  situ  sampling.  For  comparison 
purposes,  we  extracted  the  IWC  and  DE  values  from 
the  3D  satellite-inferred  cloud  fields  using  the  lati¬ 
tude  and  longitude  of  AVHRR  pixels  and  the  flight 
height  and  track.  We  then  calculated  the  mean  val¬ 
ues  for  the  pixels  associated  with  each  loop.  The  right 
panels  in  Fig.  1  illustrate  comparison  between  collo¬ 
cated  satellite-retrieved  and  in  situ  results.  The  re- 
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Fig.  1.  (Color  online)  Left  panel,  3D  IWC,  0-0.1  g/m3  and  mean  effective  ice  crystal  size  (DE,  0-200  p.m)  determined  from  a  combination 
of  the  optical  depth  and  DE  retrieval  from  the  0.63  and  3.7  pm  AVHRR  channels  aboard  the  NOAA-14  satellite  and  the  IWC  and  DE 
determined  from  the  temperature  sounding  over  the  ARM  SGP  CART  site  at  2023  UTC  on  9  March  2000.  The  3D  IWC  and  DE  results 
are  presented  in  xy,  yz,  and  xz  planes  over  a  120  km  X  40  km  X  9  km  domain.  Right  panels,  comparison  of  the  3D  DE  and  IWC  values 
inferred  from  the  satellite-temperature  and  satellite-mmCR  approaches  and  in  situ  ice  crystal  measurements  taken  onboard  the  Univer¬ 
sity  of  North  Dakota  Citation  for  the  9  March  case. 


trieved  mean  DEs  varied  between  70  and  90  pm, 
larger  than  those  derived  from  in  situ  measurements 
at  the  height  above  8  km,  but  generally  within  the 
range  of  observations  below  8  km.  Also,  the  computed 
mean  IWC  based  on  the  measured  ice  crystal  size 
distribution  varies  between  0.01  and  0.02  g/m3  and 
is  in  general  agreement  with  those  derived  from 
in  situ  measurements  in  terms  of  both  magnitude 
and  vertical  pattern.  Differences  between  the  3D 
horizontal  means  and  in  situ  values  result  from  the 
fact  that  the  latter  are  instantaneous  and  localized 
measurements.  For  comparison  purposes,  Fig.  1  also 
shows  the  mean  values  of  DE  and  IWC  as  functions  of 
AQ:  h  the  aircraft  height  inferred  from  the  satellite-mmCR 
approach. 3'lfl  The  retrieved  DE  and  IWC  vary  be¬ 
tween  60  and  100  pm  and  between  0.01  and  0.02 
g/m3,  respectively,  also  in  the  bulk  range  of  those 
derived  from  in  situ  measurements  and  in  line  with 
the  satellite-temperature  retrieval  results.  However, 
the  satellite-mmCR  approach  appears  to  be  unable  to 
produce  the  characteristics  of  increasing  DE  toward 
the  cloud  base. 

On  18  April  1997,  an  area  of  single-layer  cirrus 
cloud  was  hovering  over  the  SGP  Central  Facility  site 
at  the  NOAA-14  overpass.  This  case  is  particularly 
suitable  for  testing  of  the  3D  cirrus  mapping  scheme 
developed  above  because  the  dimension  of  this  cloudy 
area  is  about  the  size  of  a  mesoscale  model  domain, 
and  at  the  same  time,  there  were  collocated  and  co¬ 
incident  satellite  and  cloud  microphysics  data.  The 
satellite-retrieved  cloud  optical  depth  and  ice  crystal 
mean  effective  size  over  a  0.4°  X  2.0°  domain  around 
the  Central  Facility  site  reveal  the  area  of  relatively 
thick  cirrus  clouds  extends  from  about  97.2°  W  to 
98°  W.  A  small  area  of  thin  cirrus  exists  around 
96.5°  W.  The  retrieved  horizontal  optical  depth  and 
DE  range  from  0  to  8  and  from  0  to  150  pm,  respec¬ 


tively.19  The  sounding  profiles  over  the  Central  Fa¬ 
cility  at  about  2000  UTC  show  that  the  wind  direction 
was  from  the  west  and  the  speed  is  relatively  con¬ 
stant  ( — 20  m/s)  between  7  and  11  km.  Using  the 
cloud-level  temperature  sounding  and  through  the 
parameterization  described  above,  the  vertical  IWC 
and  ice  crystal  mean  size  can  be  estimated.  The 
0.4°  X  2.0°  retrieval  domain  depicted  in  Fig.  2  is  F2 
composed  of  41  scan  lines  of  AVHRR  data  in  which 
each  scan  line  contains  180  pixels.  For  the  construc¬ 
tion  of  3D  cloud  fields,  we  narrowed  the  cloud  field 
domain  to  40  km  and  selected  the  AVHRR  data  from 
a  box  of  0.4°  X  0.4°.  The  results  are  presented  in 
xy,  yz,  and  xz  planes  over  a  40  km  X  40  km  X  4.5  km 
AVHRR  data  domain,  as  shown  in  Fig.  2.  Substantial 
variability  in  IWC  and  DE  is  shown  in  the  xz  plane  as 
compared  to  those  in  the  yz  plane  associated  with  the 
prevailing  wind  direction  (parallel  to  the  y  direction), 
and  DE  is  shown  to  increase  toward  the  cloud  base. 

The  3D  IWC  and  DE  fields  for  this  case  were  com¬ 
pared  to  the  ice  crystal  size  distribution  and  IWC 
independently  derived  from  the  simultaneous  in  situ 
measurements  by  2D  probes  onboard  the  University 
of  North  Dakota  Citation.  The  flight  track  between 
2018  and  2100  UTC  consisted  of  a  racetrack  pattern 
in  six  ascending  legs  centered  at  the  ARM  central 
facility  site  near  Lamont,  Oklahoma.  This  period 
coincided  with  the  NOAA-14  AVHRR  overpass 
(—2023  UTC).  The  satellite  and  aircraft  data  extrac¬ 
tion  procedures  are  similar  to  those  for  the  9  March 
case.  The  right  panels  in  Fig.  2  illustrate  that  the 
retrieved  mean  DEs  varied  between  60  and  120  |xm, 
generally  larger  than  those  derived  from  in  situ  mea¬ 
surements  by  less  than  10  pm.  Also,  the  retrieved 
mean  IWCs  varied  between  0  and  0.3  g/m3,  in  gen¬ 
eral  agreement  with  those  derived  from  in  situ  mea¬ 
surements.  In  addition,  it  is  noted  that  the  DE  and 
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Fig.  2.  (Color  online)  Left  panel,  3D  IWC  (0-0.28  g/m3)  and  mean  effective  ice  crystal  size  (DE,  0-196  |xm)  determined  from  a  combination 
of  the  optical  depth  and  DE  retrieval  from  the  0.63  and  3.7  pm  AVHRR-NOAA  channels  and  the  vertical  IWC  and  DE  inferred  from  a 
parameterization  of  temperature  and  cloud  microphysics  over  the  ARM  SGP  CART  site  at  2023  UTC  on  18  April  1997.  The  3D  results  are 
presented  in  xy,yz,  and  xz  planes  over  a  40  km  X  40  km  X  4.5  km  domain.  Right  panels,  comparison  of  the  averaged  3D  IWC  and  DE  values 
inferred  from  the  satellite-temperature  and  satellite-mmCR  approaches  to  those  derived  from  in  situ  ice  microphysics  measurements  from 
the  University  of  North  Dakota  Citation  as  a  function  of  aircraft  leg  height  for  the  18  April  case. 


IWC  profiles  from  the  satellite-mmCR  approach 
agree  well  with  those  determined  from  the  current 
method  and  in  situ  measurements. 

In  the  preceding  comparisons,  differences  between 
the  3D  IWC  and  DE  values  inferred  from  satellite 
retrieval  and  in  situ  analysis  could  be  attributed  to  a 
number  of  approximations  made  in  the  development 
of  the  remote  sensing  algorithm  as  well  as  the  inher¬ 
ent  uncertainty  in  geographical  collocation  and  tem¬ 
poral  coincidence  between  two  data  sources.  The 
similarity  between  the  retrieved  and  in  situ  derived 
IWC  and  DE  values,  however,  is  quite  encouraging. 
Because  of  the  limitation  of  the  in  situ  measurement 
of  small  ice  crystals32  and  the  uncertainty  of  satellite 
retrieval  of  thin  clouds,33  we  anticipate  that  the  3D 
cloud  remote  sensing  technique  developed  herein  is 
most  applicable  to  cirrus  clouds  with  optical  depths 
greater  than  ~0.5  that  contain  ice  crystals  larger 
than  ~20  pm. 

3.  Three-Dimensional  Model  for  Laser 
Transmission-Backscattering  through 
Inhomogeneous  Clouds 

A.  Theoretical  Basis 

We  consider  3D  direct  transmission  and  the  trans¬ 
mission  produced  by  multiple  scattering  along  the 
laser  beam.  Figure  3(a)  depicts  a  3D  laser-target  ge¬ 
ometry  in  which  the  aircraft  is  at  the  edge  of  and 
above  an  assumed  rectangular  slab  of  a  cloud  defined 
by  a  3D  dimension  {xc,yc,zc).  The  laser-target  scan 
geometry  is  specified  by  a  vertical  scan  angle,  0, 
which  is  the  elevation  angle  of  the  laser  beam,  and  a 
horizontal  azimuthal  angle,  <J>,  which  is  the  angle 


between  the  flying  direction  of  the  aircraft  and  the 
projected  direction  of  the  laser  beam  on  the  xy  plane. 
The  positions  of  the  airborne  laser  and  the  target  are 
specified  as  (xa,ya,za)  and  (xm,ym,zj,  respectively. 
The  coordinates  of  the  incident  (xh  y„  z,)  and  exit 
(xc,  ye,  ze)  points  of  a  laser  beam  at  cloud  boundaries 
can  be  determined  from  the  laser-target  scan  geom¬ 
etry.  The  direct  transmission  through  a  3D  inho¬ 
mogeneous  cloudy  atmosphere  follows  exponential 
attenuation  and  can  be  written  in  terms  of  the 
transmitted  power  Fd  in  the  form 


Fd  =  F0  exp 


-t* 


pe(s')ds' 


(3.1) 


where 


U  Pnc.a^  "f  Pnc,6^U  Pe  Pair  Pacr  T  Pdu . 

In  the  above  equations,  F0  is  the  laser  power  in  watts; 
the  total  extinction  coefficient  p,,  is  the  sum  of  the 
scattering  coefficient  of  air  molecules  pair,  the  extinc¬ 
tion  coefficient  of  aerosols  paor,  the  absorption  coeffi¬ 
cient  of  water  vapor  times  water  vapor  density  kup, 
and  the  extinction  coefficient  of  cloud  particles  pcU 
along  the  path  of  the  laser  beam;  and  the  exponential 
attenuation  outside  the  cloud  is  also  accounted  for  in 
Eq.  (3.1),  where  pnc/J  and  pnci  are  extinction  coeffi¬ 
cients  due  to  noncloud  species  (air  molecules,  aero¬ 
sols,  and  water  vapor)  above  and  below  the  cloud, 
respectively.  The  parameter  u  denotes  the  distance 
BTg  between  the  cloud  base  and  the  target  and  v 
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Fig.  3.  (Color  online)  (a)  Laser  transmission  through  a  3D  cloud. 
The  positions  (x„,  y„ ,  z„),  (x,„,  y,„,  z,„),  (xh  yh  z;),  and  fo,  ye,  zr)  denote 
coordinates  of  the  airborne  laser,  target,  incident  point,  and  exit 
point;  yc,  z,.)  is  the  3D  cloud  dimension,  where  zr  =  z,  —  z*;  0  is 
a  laser  scan  angle;  <}>  is  a  horizontal  azimuthal  angle; is  an  angle 
subtended  by  the  target  at  point  E;  r„,  is  the  effective  radius  of  the 
target;  s'  =  EB;  u  =  BTg ;  v  =  AT;  and  s  -  s'  =  ET.  (b)  Laser 
backscattering  through  a  3D  cloud.  'I/,  is  the  angle  subtended  by 
the  aircraft  at  point  E,  ra  is  the  effective  radius  of  the  aircraft, 
s  -  sh'  =  ETg ,  and  sh'  =  EA.  (c)  Laser  transmission  through  a  3D 
inhomogeneous  cloud  in  a  spherical  atmosphere.  All  symbols  are 
the  same  as  in  (a),  except  that  s*  =  TlT2,  s'  =  ET2,  u  =  T2Tg,  and 
'Ll  5'  is  an  angle  subtended  by  the  target  at  point  E. 


denotes  the  distance  AT  between  the  aircraft  and  the 
cloud  top.  In  the  present  study,  (3aer  is  obtained  from 
the  MODTRAN  default  profile,  and  p  is  taken  from 
the  U.S.  Standard  Atmosphere. 

In  reference  to  Fig.  3(a),  the  parameter  s'  is  the 
distance  along  the  laser  beam  between  a  point  ( x ,  y,  z) 
in  the  cloud  and  the  exit  point  at  the  cloud  boundary 
defined  by 


s'  =  (z-zb)/ p,,  (3.2) 

where  p,  is  sin  1 0  | ,  zh  is  the  cloud-base  height,  and  0 
is  the  scan  angle  of  the  laser  beam  in  the  xz  plane 
given  by 


<t>  =  tan  l[(ym-ya)/{xm-xa)\  (3.4) 

The  distance  along  the  laser  beam  between  the  cloud 
top  and  base  is  given  by 

s  =  (zt-  zb)/|x  =  {(xc  -  x,.)2  +  (y,  -X)2  +  (ze  -  Zi)2}1'2, 

(3.5) 

where  zt  is  the  height  of  the  cloud  top.  The  points 
(Xi>  y <n  Ai)*  (Ai  y mi  ** m)i  (-A  y i i  a)  and  (xe,  ye,  zr)  are  the 
coordinates  of  the  aircraft,  the  target,  and  the  inci¬ 
dent  and  emergent  points  of  the  laser  beam  on 
the  cloud,  respectively.  The  parameters  xc,  yc,  and  zc 
denote  the  3D  cloud  dimensions. 

Formulation  of  the  first-order  scattering  in  the  3D 
case  is  similar  to  the  2D  case  and  can  be  written  in 
the  form 


F(1)(0,  fit)  =  exp(— tc)  j  J°\s',  CL) 
J  0 


X  exp 


Pt,(s")ds" 


Pc(s')ds'. 


(3.6) 


However,  the  path  lengths  s’  and  s  are  functions  of 
3D  cloud  geometry  defined  in  Eqs.  (3.2)— (3.5).  Based 
on  the  successive-order-of-scattering  approach,  the 
first-order  source  function  for  forward  scattering  is 
given  by 


m(s'). 


J°\s',  Cl)  =  -^F<°>(s\  Cl)  j  Pcld(0,  s') sin  0d0, 

(3.7a) 


Vi 


where  the  term 


F(0)(s',  Cl)=F0exp 


Pc(s")ds" 


(3.7b) 


The  single-scattering  albedo,  Trr(s'),  and  the  phase 
function  of  cloud  particles,  Pckl,  are  functions  of  s' 
involving  3D  geometry,  and  the  angle  T,  within 
which  scattering  occurs  is  given  by  tan_1[r„,/(« 
+  s')]-  Note  that  the  phase  function  for  aerosols  (and 
air  molecules)  in  forward  directions  is  more  than  two 
orders  smaller  than  that  for  ice  particles,  and  for 
practical  purposes  the  former  can  be  neglected  in 
forward-scattering  calculations. 

In  a  manner  similar  to  the  first-order  forward¬ 
scattering  formulation,  the  source  function  that  is 
due  to  the  first-order  backscattering  can  be  expressed 
in  the  form 


0  =  tan  '[cos  -  z„)/ ( xm  -  x„)],  (3.3) 

with  4>  being  the  azimuthal  angle  of  the  laser  beam: 


J°\sb',  O)  =  ^F<0>(s6',  Cl) 
X  sin  0d0, 


-n 


PM®,  V) 
(3.8a) 
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where  the  term 


craft  and  target  are  above  the  cloud,  F,  is  defined  by 


Fm(sb',  fl)  =F0  exp 


pc(S")ds" 


0 


(3.8b) 


Fs  =  exp 


pe(s')ds' 


0 


(3.12) 


st,'  is  the  distance  from  the  aircraft  to  a  point  along 
the  laser  beam,  as  shown  in  Fig.  3(b);  and  Pavc(0,  sbj 
is  the  phase  function  that  accounts  for  both  cloud  and 
aerosol  particles.  For  backscattering,  contribution 
from  aerosols  becomes  significant  when  thin  cirrus 
clouds  are  involved,  and  the  angle  4ft  is  given  by 
tan ~l{ra/sbj,  where  r„  denotes  the  effective  radius  of 
the  aircraft.  Note  that  the  phase  functions  for  aero¬ 
sols  and  air  molecules  are  of  the  same  order  as  that 
for  cloud  particles  and  must  be  accounted  for  in  the 
analysis.  Thus  formulations  for  the  first-order  for¬ 
ward  [Eqs.  (3.6),  (3.7a),  and  (3.7b)]  and  backward 
scattering  [Eqs.  (3.8a)  and  (3.8b)]  differ  not  only  in 
geometry  but  also  in  the  use  of  the  phase  function. 

Furthermore,  using  the  principle  of  successive  or¬ 
der  of  scattering  along  the  laser  beam,  the  nth-order 
scattering  can  be  expressed  by 


F(,,)( 0,  fl)  =  exp(-Tc)  fl) 


x  exp 


pc,(s")ds" 


3t,(s')ds',  (3.9) 


where  the  source  function 


TiT(S  ) 

fi)  =  I  V,  fi)Pdd(0,  s') 

X  sin  0d0,  (3.10) 

where  0„  denotes  the  scattering  angle  within  which 
the  nth-order  scattering  contributes  to  the  source 
function  integral,  and  4ft  is  an  angle  associated  with 
0„  determined  by  the  scattering  geometry.  Likewise, 
the  nth-order  backscattering  can  also  be  expressed  in 
a  similar  form. 

The  distance  between  the  aircraft  and  the  engaging 
target  is  generally  of  the  order  of  hundreds  of  kilo¬ 
meters.  To  accurately  simulate  laser  transmission 
along  such  a  long  path  length,  the  sphericity  of  the 
cloud-top  surface  due  to  the  Earth’s  curvature  must 
be  taken  into  account.  As  illustrated  in  Fig.  3(c),  we 
consider  the  geometry  where  both  the  aircraft  and 
the  target  are  above  the  cloud.  The  direct  transmis¬ 
sion  of  a  laser  beam  in  a  3D  spherical  geometry,  Fd*, 
is  similar  to  Eq.  (3.1)  and  can  be  expressed  by 

Fd*  =  FdFs,  (3.11) 

where  Fd  is  the  direct  transmission  for  the  clear-sky 
line  of  sight,  and  Fs  is  an  additional  cloud  transmission 
produced  by  spherical  geometry.  When  both  the  air- 


where  s*  is  the  distance  along  the  laser  beam  be¬ 
tween  the  incident  (x„  yh  zj  and  exit  (xe,  ye,  zj  points 
at  the  cloud  top,  and  s'  is  the  distance  along  the  laser 
beam  between  a  point  in  the  cloud  and  the  exit 
(xc,  y„  zj  point  at  the  cloud  top.  Forward  scattering  is 
the  sum  of  the  forward  scattering  produced  in  the 
plane-parallel  geometry,  Fm(0,  fl),  plus  the  forward 
scattering  in  the  portion  of  the  cloudy  atmosphere 
that  is  related  to  the  spherical  geometry: 

F*O)(0,  fl)  =F(1)(0,  fl)  +FS(1)(0,  fl).  (3.13) 

On  the  basis  of  the  successive  order-of-scattering  ap¬ 
proach,  the  second  term  is  given  by 


Psfl>(0,  f!)  =  exp(-Tc)  fl) 


X  exp 


pc(s'')ds" 


pt.(s')ds',  (3.14) 


where  s*  is  the  distance  between  the  incident  and 
exit  points.  The  expression  for  the  source  function 
J°V,  fl)  in  Eq.  (3.14)  is  the  same  as  that  defined  in 
Eqs.  (3.7a)  and  (3.7b),  except  that  the  angular  inte¬ 
gration  limits  are  determined  by  a  new  angle  4ft'*  and 
the  path-length  integration  for  the  direct  transmis¬ 
sion  is  from  0  to  (s*  -  s').  Likewise,  the  source  func¬ 
tion  due  to  the  first-order  backscattering  can  be 
obtained  from  the  angular  and  path-length  integra¬ 
tion.  Furthermore,  the  nth-order  scattering  terms  for 
transmission  and  backscattering  are  the  same  as 
those  defined  in  Eqs.  (3.9)  and  (3.10),  except  for  the 
angular  and  path-length  integration  limits  that  fol¬ 
low  the  spherical  geometry  formulation. 

B.  Application  of  the  Satellite-Inferred  Three-Dimensional 
Results  to  Three-Dimensional  Laser  Transmission  and 
Backscattering 

We  investigate  the  sensitivity  of  laser  transmission 
and  backscattering  on  variation  in  the  cloud  optical 
depth,  particle  size,  and  cloud-base  height,  as  well  as 
the  positions  of  the  aircraft  and  target.  In  particular, 
we  study  deviation  of  the  2D  transmission  and  back- 
scattering  results  from  the  present  3D  values  using 
the  same  cloud  properties.  Here  2D  implies  that  a 
cloud  field  is  averaged  over  the  y  direction  for  laser 
transmission  and  backscattering  formulations.18 

Figure  4  shows  three  vertical  transects  aty  =  1,  20, 
and  39  km  corresponding  to  the  3D  cirrus  cloud  IWC 
and  DE  determined  for  the  18  April  1997  case  dis¬ 
cussed  above.  The  aircraft  (laser  beam)  is  shown  at 
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Fig.  4.  (Color  online)  3D  cirrus  IWC  and  DE  fields  mapped  aty  =  1,  20,  and  39  km  obtained  for  the  18  April  1997  case  over  the  ARM 
SGP  Central  Facility.  The  dark  areas  are  areas  with  no  observations.  Also  shown  are  the  positions  of  the  aircraft  and  the  moving  target. 


the  left  upper  edge  (x  =  0,  z  =  za)  of  these  diagrams, 
while  the  target  is  located  at  the  right  edge  (x  =  xm). 
These  transects  of  cirrus  cloud  properties  differ  signif¬ 
icantly  and  therefore  cannot  be  assumed  to  be  homo¬ 
geneous  in  the  y  direction.  For  example,  the  IWC 
transect  at  y  =  1  km  shows  an  area  of  maximum 
values  near  x  =  15  km  and  2  =  9  km,  whereas  at 
y  =  20  and  39  km,  there  is  no  such  maximum.  More¬ 
over,  the  Dc  transect  shows  a  maximum  at  y  = 
20  km. 

In  Fig.  5,  we  present  the  mean  and  standard  devi¬ 
ation  of  3D  results  for  direct  transmission  and  first- 
order  forward  and  backward  scattering  based  on  the 
five  vertical  transects  (1,  11,  20,  29,  and  39  km)  as 
functions  of  target  height.  The  standard  deviation  is 
denoted  by  a  horizontal  bar.  For  comparison  pur¬ 
poses,  the  2D  results  are  also  displayed  (dotted 
curves).  Substantial  differences  are  shown  between 
2D  and  the  3D  mean,  particularly  in  the  cloud  region 


for  larger  optical  depth.  For  direct  transmission  and 
first-order  scattering,  standard  deviations  are  larger 
inside  the  cloud  since  the  path  length  within  the 
cloud  is  longer  and  the  laser  transmission  is  strongly 
subject  to  cloud  inhomogeneity.  For  first-order  back- 
scattering,  standard  deviations  in  the  cloud  are 
smaller,  since  it  is  primarily  determined  by  the  cloud- 
top  portion  near  the  aircraft.  The  total  transmission, 
the  sum  of  direct  transmission  and  first-order  scat¬ 
tering,  may  be  expressed  in  a  parameterized  form, 
exp(-T*/p-),  where  t*(~0. 95-0.99)  is  an  adjusted 
optical  depth  accounting  for  both  contributions. 
Backscattering  depends  on  the  characteristics  of  the 
particles’  phase  function  and  the  intricate  geometry 
involving  laser  and  target.  For  this  reason,  develop¬ 
ment  of  an  analytical  model  for  physical  explanation 
is  not  straightforward. 

Finally,  the  effects  of  the  Earth’s  sphericity  on  laser 
transmission  and  backscattering  were  evaluated.  Fig- 
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3D  Average  2D  Inhomgeneous  Cloud 


Forward  scattering  (W) 


10'  103  to '  10'  10s  10’  10'  10’  10'  10’  10°  10' 

Backscattering  (W) 

Fig.  5.  (Color  online)  Direct  transmission,  forward  scattering, 
and  backscattering  for  a  2D  inhomogeneous  (averaged  over  the  y 
direction)  and  moan  values  of  those  for  the  3D  cirrus  cloud  field 
based  on  the  five  vertical  transects  (1,  11,  20,  29,  and  39  km)  are 
displayed.  Standard  deviations  of  the  3D  mean  values  based  on  the 
five  3D  transect  values  are  also  plotted  in  terms  of  horizontal  bars. 
X  =  1.315  p.m,  F0  =  106  W,  RK  =  0.2,  r„  =  4  m,  z0  =  11.2  km, 
zt  =  11.1  km,  and  z*  =  7.5  km. 


ure  6  shows  plots  for  the  3D  direct  transmission  and 
first-order  forward  scattering  and  backscattering 
through  100  km  X  100  km  X  4.5  km  and  200  km 
X  200  km  X  4.5  km  cloud  fields  embedded  in  spher¬ 
ical  and  plane-parallel  atmospheres.  The  3D  cloud 
horizontal  fields  of  40  km  X  40  km  are  adjusted  uni¬ 
formly  to  fit  in  the  100  km  X  100  km  and  200  km 
X  200  km  dimensions  for  the  purpose  of  understand¬ 
ing  the  effect  of  the  Earth’s  curvature  on  laser  trans¬ 
mission.  For  direct  transmission,  the  effect  of  a 
spherical  Earth  is  more  significant  for  a  longer  hori¬ 
zontal  distance  of  200  km  than  100  km,  particularly 
when  the  target  height  is  between  the  cloud  top  and 
2  km  above  it,  where  transmission  can  be  signifi¬ 
cantly  reduced  (by  two  orders  of  magnitude)  by  the 
Earth’s  curvature.  The  forward  scattering  is  en¬ 
hanced  near  the  cloud  top  due  to  the  additional  scat¬ 
tering  by  laser  transmission  through  the  curved 
cloud-top  portion.  The  slight  fluctuation  above  the 


3D  Spherical  3D  Plane  Parallel 


Direct  Transmission  (W) 


Forward  scattering  (W) 


10‘  103  10'  10'  10’  10’  10“*  10 3  10'  10’  10”  10' 


Backscattering  (W) 

Fig.  6.  (Color  online)  Comparisons  of  the  direct  laser  transmis¬ 
sion,  forward-scattering  transmission,  and  backscattering  for 
cloud  fields  in  plane-parallel  and  spherical  atmospheres  assuming 
horizontal  distances  of  100  and  200  km.  The  3D  cloud  fields 
(40  km  x  40  km  in  the  horizontal)  are  adjusted  to  fit  in  the 
100  km  X  100  km  and  200  km  X  200  km  dimensions. 


cloud  is  due  to  localized  strong  scattering  near  the 
cloud  top.  The  effects  of  the  spherical  geometry  on  the 
first-order  backscattering  and  reflection  from  the  tar¬ 
get  are  small  because  of  small  backscattering  values. 
For  both  Figs.  5  and  6,  the  results  do  not  show  sen¬ 
sitivity  to  cloud  inhomogeneity  above  z  =  12  km. 
However,  we  have  included  the  curves  above  this 
height  to  be  consistent  with  the  format  in  our  previ¬ 
ous  papers  published  in  this  journal.17-38 

4.  Concluding  Remarks 

We  have  developed  an  approach  for  the  determination 
of  3D  IWC  and  DE  fields  of  cirrus  clouds  based  on  a 
combination  of  satellite-retrieved  optical  depth  and  DE 
in  a  horizontal  plane  and  the  vertical  IWC  and  DE 
determined  through  a  parameterization  of  cloud  mi¬ 
crophysics  and  temperature.  To  investigate  the  reli¬ 
ability  of  this  parameterization  approach  for  the 
mapping  of  3D  cloud  fields,  we  have  selected  AVHRR- 
NOAA  data  associated  with  the  single-layer  cirrus 
clouds  that  occurred  at  the  ARM  SGP  CART  site  in  the 
northern  Oklahoma  area  on  18  April  1997  and  9 
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March  2000  during  a  Spring  Cloud  IOP.  Retrieval  and 
analysis  were  carried  out  to  construct  3D  IWC  and  DE 
fields  from  the  AVHRR  data,  which  were  then  com¬ 
pared  to  collocated  and  coincident  in  situ  measure¬ 
ments  of  the  ice  crystal  size  distribution  and  IWC 
independently  derived  from  the  optical  probes  onboard 
the  University  of  North  Dakota  Citation.  The  satellite- 
retrieved  3D  IWC  and  DE  values  show  a  general  and 
consistent  agreement  with  the  concurrent  and  collo¬ 
cated  airborne  microphysics  measurements.  The  accu¬ 
racy  of  small  ice  crystals  in  both  in  situ  measurements 
and  satellite  remote  sensing  is  a  subject  requiring  fur¬ 
ther  investigation. 

A  3D  laser  transmission  and  backscattering  model 
for  inhomogeneous  clouds,  particularly  applicable  to 
high  thin  cirrus,  has  been  developed  for  an  aircraft- 
target  system.  We  approached  this  problem  by  using 
the  principle  of  the  successive  order  of  scattering  for 
radiative  transfer  in  which  the  attenuation  and  scat¬ 
tering  of  a  laser  beam  can  be  traced  in  accordance  with 
the  scattering  geometry  for  clouds  in  plane-parallel 
and  spherical  atmospheres.  For  transmission,  the  di¬ 
rect  component,  the  first-order  scattering,  and  the  nth- 
order  multiple  scattering  were  formulated  by  means  of 
the  source-function  technique  in  the  context  of  the 
successive-order-of-scattering  approach.  Similar  for¬ 
mulations  were  carried  out  for  backscattering  involv¬ 
ing  the  laser  energy  available  at  the  position  of  the 
aircraft.  We  performed  numerical  simulation  of  the 
transmitted  and  backscattered  laser  energies  for  a  3D 
inhomogeneous  cirrus  cloud  field.  The  simulation  re¬ 
sults  for  3D  direct  transmission,  forward-scattered 
transmission,  and  backscattering  have  been  compared 
to  those  for  the  2D  model  and  substantial  differences 
are  shown.  The  effect  of  the  Earth’s  curvature  becomes 
particularly  important  for  long  path  transmission  in¬ 
volving  the  airborne  laser  and  moving  target  engage¬ 
ment. 

In  the  present  study,  we  have  used  the  AVHRR- 
NOAA  data  for  the  construction  of  3D  cloud  fields  in 
view  of  the  availability  of  collocated  and  coincident 
airborne  ice  microphysics  measurements  for  validation 
purposes.  However,  the  parameterization  approach 
developed  in  this  paper  can  be  directly  applied  to  the 
data  sets  from  MODIS  onboard  the  NASA  Terra  and 
Aqua  satellites  as  well  as  the  future  visible-infrared 
AQ:  i  imager  radiometer  suite  channels  to  be  included  in  the 
National  Polar-orbiting  Operational  Environmental 
Satellite  System  platform  for  the  3D  cirrus  cloud 
study.  The  availability  of  these  satellite  data  for  3D 
cloud  mapping  can  be  directly  employed  in  connection 
with  the  simulation  of  3D  airborne  laser  transmission 
and  backscattering  in  spherical  atmospheres  for  tacti¬ 
cal  application. 
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